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Non-Markovian effects on the dynamics of entanglement 
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A procedure that allows to obtain the dynamics of A'^ independent bodies each locally interacting 
with its own reservoir is presented. It relies on the knowledge of single body dynamics and it is 
valid for any form of environment noise. It is then applied to the study of non-Markovian dynamics 
of two independent qubits, each locally interacting with a zero temperature reservoir. It is shown 
that, although no interaction is present or mediated between the qubits, there is a revival of their 
entanglement, after a finite period of time of its complete disappearance. 
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Entanglement is relevant to different fundamental [l|,0| 
aspects of quantum theory andpractical aspects of quan- 
tum information processing [3|. Recently much inter- 
est has arisen in the evolution of the joint entanglement 
of a pair of qubits exposed to local noisy environments. 
The reason is related to the discovery by Yu and Eberly 
Q that for this system, rather surprisingly, the Marko- 
vian dynamics of the joint qubits entanglement and single 
qubit decoherence may be rather different. The aspect 
that has mostly drawn attention is the possibility of a 
complete disappearance of entanglement at finite times. 
The occurrence of this phenomenon, termed "entangle- 
ment sudden death" (ESD), has been shown in a quan- 
tum optics experiment Q . The intrinsic interest and po- 
tential importance of ESD, for example in the application 
range of quantum error correction methods, has lead to a 
flow of analysis that stud y it s appearance under different 
circumstances 0, 0, H, i,M [ill 111, [H ■ 

Disentanglement is related to the birth of body- 
environment correlations. It is therefore of interest to 
investigate the role played on its evolution by non- 
Markovian effects. In fact, although Markovian dynamics 
includes a level of back-reaction, it neglects the entangle- 
ment that arises between bodies and bath modes during 
the evolution. Alt houg h some work has treated of this 
aspect 14, [l^ 16, 17, it should be considered an 
attractive theoretical challenge to extend the results ob- 
tained under various conditions in the Markovian regime 
to the non-Markovian case [ll]. The aim of this letter 
is to address this point first by adopting a procedure to 
obtain the dynamics of N independent bodies, locally in- 
teracting with reservoirs and without restriction on the 
nature of environmental noise, if the single body dynam- 
ics is known. We shall then use this approach to explicitly 
investigate the entanglement dynamics of two qubits lo- 
cally interacting with a zero temperature non-Markovian 
environment. 

To describe the method we consider a system formed 
by two non-interacting parts A,B, each part consisting 
of a qubit S — A,B locally interacting respectively with 
a reservoir Rs = Ra,Rb- Each qubit and the corre- 



sponding reservoir are initially considered independent. 
For each part, the reduced density matrix evolution for 
the single qubit S — A,B is given by 

= Tr;^„ {uHt)pHO) ^ p''HO)U^Ht)} , (1) 

where the trace is over the reservoir Rs degrees of free- 
dom and U^{t) is the time evolution operator for the part 
S. In terms of the Kraus operators p{t), the former 
equation becomes [2^ 



(2) 



The assumption of independent parts implies that the 
time evolution operator U'^{t) of the complete system 
f = A + B factorizes as U^{t) = U^{t) ® U^{t). It fol- 
lows that the Kraus representation of the reduced density 
matrix for the two-qubit system T = A -\- B reads like 

(3) 

Given the basis {|0), |1)} for each qubit, inserting unity 
operators / = |0) (0| -I- 11) (1| between Kraus operators and 
density matrices in Eq. it follows that the dynamics 
of each qubit has the form 

W mm' 

(4) 

Adopting the same procedure for p^ in the form of 
Eq. ([3]), the dynamics of the two-qubit system is given 

by 



Pii' (0 



^ 4''(Oi3,T'(Opn',w(o), (5) 



where the indexes i,j,l,m — 0,1. Eqs. (g]), (O clearly 
show that the dynamics of two-qubit density matrix el- 
ements can be obtained by knowing that of the single 
qubit. The validity of above procedure can be straight- 
forwardly extended to any multipartite and multilevel 
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system (qudit), provided that the different parts, "qudit 
+ reservoir", are independent. 

We now apply the results obtained above to study non- 
Markovian effects on the entanglement dynamics of two 
qubits, each interacting only and independently with its 
local environment. To this aim we shall consider the 
single "qubit+reservoir" Hamiltonian given by 



k 



(6) 



with B — J2k9kbk, where ojo is the transition frequency 
of the two-level system (qubit) and a± are the system 
raising and lowering operators while the index k labels 
the field modes of the reservoir with frequencies uj^, bj_, 
are the modes creation and annihilation operators and gk 
the coupling constants. The Hamiltonian of Eq. ^ may 
describe various systems as for example a qubit formed by 
an exciton in a potential well environment. However to 
fix our ideas we shall take it to represent a qubit formed 
by the excited and ground electronic state of a two-level 
atom interacting with the reservoir formed by the quan- 
tized modes of a high-Q cavity. At zero temperature, this 
Hamiltonian represents one of the few open quantum sys- 
tems amenable to an exact solution [2l| . The dynamics of 
qubit S is known to be described by the reduced density 
matrix E^E^ 



pmpt 



p'{t) =1 I , (7) 



where the function Pt obeys the differential equation 



Pf 



dhf{t - h)Pt, , 



(8) 



and the correlation function f{t — ti) is related to the 
spectral density J{i^) of the reservoir by 



f{t — ti)— / dujJ{uj)exp[i{ujo — iij){t — ti 



(9) 



The exact form of Pt thus depends on the particular 
choice for the spectral density of the reservoir. Because 
the Hamiltonian of Eq. ([6]) represents a model for the 
damping of an atom in a cavity, we consider then the case 
of a single excitation in the atom-cavity system. For the 
effective spectral density J(w), we take the spectral dis- 
tribution of an electromagnetic field inside an imperfect 
cavity supporting the mode loq, resulting from the com- 
bination of the reservoir spectrum and the system reser- 
voir coupling with 70 related to the microscopic system- 
reservoir coupling constant, of the form [2^ 



1 



7oA2 



2tt (loo - cj)2 + A2 



(10) 



The correlation function ^ corresponding to the spec- 
tral density of Eq. (jlOp has an exponential form with A as 



decay rate. The parameter A, defining the spectral width 
of the coupling, is then connected to the reservoir corre- 
lation time tb by the relation tb ~ A^^. On the other 
hand the parameter 70 can be shown to be related to the 
decay of the excited state of the atom in the Markovian 
limit of flat spectrum. The relaxation time scale t^^ over 
which the state of the system changes is then related to 
70 by TR ^ 7o"^ 

Using the spectral density of Eq. (fTO|) in the correlation 
function of Eq. in the subsequent analysis of the 
function Pt of Eq. ([S]), typically a weak and a strong 
coupling regime can be distinguished. For weak regime 
we mean the case 70 < A/2, that is tu > 2tb- In this 
regime the relaxation time is greater than the reservoir 
correlation time and the behaviour of Pt is essentially 
a Markovian exponential decay controlled by 79. In the 
strong coupling regime, that is for 70 > A/2, or tr < 2tb, 
the reservoir correlation time is greater or of the same 
order of the relaxation time and non-Markovian effects 
become relevant. For this reason we are interested in 
this regime and we shall mainly limit our considerations 
to this case. Within this regime, the function Pt assumes 
the form 
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(11) 



where d = \/27oA — A^. Pt presents oscillations describ- 
ing the fact that the decay of the atom excited state is in- 
duced by the coherent processes between the system and 
the reservoir. In particular he function Pt has discrete ze- 
ros at i„ = 2 [riTT — arctan(rf/A)] /d, with n integer. We 
note that the solution in the weak coupling regime can be 
obtained by the former one simply substituting the har- 
monic functions with the corresponding hyperbolic ones 
and d with id. 

Now we are ready to use, following the procedure de- 
scribed before, the evolution of the reduced density ma- 
trix elements for the single qubit to construct the reduced 
density matrix jiF for the two-qubit system. In the stan- 
dard product basis B = {|1) = |11),|2) = |10),|3) = 
|01), |4) EE |00)}, using Eqs. (H), © and dH), we obtain 
the diagonal elements 



Pl2{t) 

pUt) 
pLit) 



- P22(0)Pt+pil(0)Pt(l 

= pUo)Pt + pl,Pt{i-Pt). 

= l-bn(0+P^2(t)+pl3(i)], 



(12) 



and the non-diagonal elements 



Pl2{t) 

pIS) 
pUt) 

pUi) 



= plMPr\ pUt)-plMPf'\ 
^ plMPt, pLit) = p^3{o)Pt, 

= ^/Pt{p^M+p^3(.m-Pt)], 



(13) 
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and pfj{t) = pj*{t), where p'^(t) is a hermitian matrix. 
In order to follow the entanglement dynamics of the bi- 
partite system, we use Wootters concurrence jU- This 
is obtained from the density matrix fp" for qubits A and 
B as CpT (t) = max{0, y/Xi — ^/X2 — a/As — \A4}, where 
the quantities Aj are the eigenvalues of the matrix C 

C = p^(aX)p^*(aX), (14) 

arranged in decreasing order. Here fP"* denotes the com- 
plex conjugation of in the standard basis and ay is 
the well-known Pauli matrix expressed in the same ba- 
sis. The concurrence varies from C = for a disentangled 
state to C = 1 for a maximally entangled state. 

The form of Eqs. (|12p and is such that one can 
study the entanglement evolution for any initial state. 
We shall however restrict our analysis to the initial en- 
tangled states 

1$) = a|01) -^/3|10), l^-) = a|00) +/3|11), (15) 

where a is real, f3 = \f3\e^^ and + = 1. For these 
two entangled states, the initial total density matrix has 
an "X" structure [l9| which is maintained during the 
evolution, as easily seen from Eqs. and (fTH]) . In par- 
ticular the concurrence, for the initial states of Eq. 
is given by 

C*(t) = max{0,2|p^3(0| - 2^ pf,{t)pj,{t)}, 

C^{t) = mB.^{Q,2\pUt)\-'i^pUt)pUt)}^ (16) 
and using Eqs. (fT2|) and (fTS]) we obtain 

C^{t) = max{0,2yr~^Pta}, 

C^{t) = max{0,2v/l - a2Ff[a- 1 - a"^ [l - Pt)]} . 

(17) 

The time behavior of the concurrences C$ and as 
a function of and the dimensionless quantity 7oi are 
plotted for Pt given by Eq. (fTTjl in Figs. 1 and 2 (for 
A/70 = 0.1). This is a condition that can be well within 
the current experimental capabilities. In fact, cavity 
QED experimental configurations have been realized us- 
ing Rydberg atoms with lifetimes Tat ~ 30ms, inside 
Fabry-Perot cavities with quality factors Q w 4.2 x 10^° 
corresponding to cavity lifetimes Tcav ~ 130ms |25[; these 
values correspond to 2 A/70 ~ 0.2 

Fig. 1 shows that, in the non-Markovian regime, the 
concurrence C$ periodically vanishes according to the 
zeros of the function Pt, with a damping of its revival 
amplitude. This behavior is evidently different from 
Markovian, where in contrast C$ decays exponentially 
and vanishes only asymptotically p] . The Markovian de- 
cay rate is however larger than the initial non-Markovian 
one, as shown in Fig. 3 for the maximally entangled case 
a2 = 1/2. 




FIG. 1: Concurrence for the initial state q|01) + /?|10) as a 
function of the dimensionless quantity 7ot and , in a real- 
istic CQED condition (A = O.I70). 




FIG. 2: Concurrence for the initial state q|00) + /9ill) as a 
function of the dimensionless quantity 7ot and , in a real- 
istic CQED condition (A = O.I70). 



Fig. 2 shows that the entanglement represented by Ciif 
has a similar behavior to C$ for > 1/2. In contrast for 
< 1/2 two ranges of parameter may be distinguished. 
In one there is ESD because vanishes permanently 
after a finite time, similar to the Markovian case [1, [J. 
In the second, revival of entanglement appears after pe- 
riods of times when disentanglement is complete. This 
behavior is more evident in the plot of Fig. 4 obtained 
under stronger non-Markovian conditions. This revival 
phenomenon is induced by the memory effects of the 



c(t) 




FIG. 3: Concurrence for the initial state (|01) + |10))/\/2 as a 
function of the dimensionless quantity 7ot in non-Markovian 
regime (solid line; A — O.I70) and Markovian regime (dashed 
line; A = 670). 



4 



0.5 




Yot 



FIG. 4: Concurrence for the initial state q:|00) + as a 

function of the dimensionless quantity jot for — 1/3, in 
strong couphng (A = O.OI70). 



reservoirs, which allows to the two-qubit entanglement 
to reappear after a dark period of time, during which the 
concurrence is zero. This phenomenon of revival of en- 
tanglement after finite periods of " entanglement death" 
appears to be linked to the environment-single qubit non- 
Markovian dynamics. In this sense this result differs from 
the revival of entanglement previously obtained in the 
presence of interaction among qubits or because of their 
interaction with a common reservoir 

0,1 



11, US 



The 

physical conditions examined here are moreover more 
similar to those typically considered in quantum com- 
putation, where qubits are taken to be independent and 
where qubits interact with non-Markovian environments 
typical of solid state micro devices [2^ . 

The above analysis can be easily extended to study en- 
tanglement dynamics starting from different initial condi- 
tions and to take into account finite temperature effects. 
In particular, starting from a Werner state [27| one gets 
an entanglement behavior structurally similar to that ob- 
tained in this letter for the states of Eq. (fT5)) with ESD 
and entanglement revival periods. The details of the evo- 
lution for this case and finite temperature effects for non- 
Markovian dynamics will be considered elsewhere. 

In conclusion we have presented a procedure that al- 
lows in principle to obtain the dynamics of a system of 
N independent bodies, each locally interacting with an 
environment, as long as the single system dynamics is 
known. This procedure is valid for any form of single 
body-environment interaction. It has been applied to 
the case of two qubits interacting with the environment 
where non-Markovian effects are present. In particular 
the model described has been identified with a system 
made by two atoms each in a high-Q cavity. For this 
case the Hamiltonian dynamics of the single qubit can 
be solved exactly and no problem about map positivity 
arise. We have found that non-Markovian effects influ- 
ence the entanglement dynamics and may give rise to a 
revival of entanglement even after complete disentangle- 
ment has been present for finite time periods. This ef- 



fect, arising for completely independent systems, is only a 
consequence of the non-Markovian behavior of the single 
qubit-reservoir dynamics. 

These results show that entanglement dynamics may 
present facets than one may not simply expect from single 
qubit dynamics and may also lead to the possibility of 
recovering the entanglement initially present. 
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